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Abstract 

We construct an analytic solution of the Einstein-5'C/(2)-Yang-Mills system as the 
holographic dual of an anisotropic superfluid near its critical point, up to leading correc- 
tions in both the inverse Yang-Mills coupling and a symmetry breaking order parameter. 
We have also calculated the ratio of shear viscosity to entropy density in this background, 
and shown that the universality of this ratio is lost in the broken symmetry direction. 
The ratio displays a scaling behavior near the critical point with critical exponent f3 = 1, 
at the leading order in the double expansion. 
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1 Introduction 



Gauge/gravity duality has led to many useful insights into strongly coupled field theories. 
Recently, fluid/gravity duality has been widely studied, providing useful information about 
conformal fluid dynamics, the effectively long wavelength description of conformal field theory. 
Many crucial quantities characterizing conformal fluids can be obtained as real-time Green's 
functions in dual gravitational theories jH [3] . 

The most celebrated example from holographic fluid dynamics is the ratio of shear viscosity 
n to entropy density s, 2 = [UE1E1E]- The ratio seems to be universal for theories 
with weakly coupled gravity duals E • However, in recent studies of anisotropic conformal 
fluids [TH [13], the universality of the ratio turns out to be violated [TU [17]. These studies 
have considered an exact solution of the bulk Einstein-5£/ (2)-Yang-Mills action @ for an AdS 
black brane in 5 dimensions with nonzero chemical potential. Corresponding to the chemical 
potential, the temporal component of the gauge potential is proportional to a 3 of the SU{2) 
gauge group. The boundary metric enjoys 50(3) global symmetry (symmetry under spatial 
rotations). At high temperature (equivalently small chemical potential fi), the system stays in 
the isometric phase. However at a certain chemical potential fi = fi c , the 50(3) symmetry is 
spontaneously broken to 50(2) because one of the spatial components of the Yang-Mills field 
develops a non-trivial zero mode (as a solution of the linearized Yang-Mills field equations) 
and it is normalizable solution. It turns out that in the region of fi > fj, c , this mode condenses 
and there is a new anisotropic superfluid phase. 

As long as we are looking at a solution with 50(3) symmetry, the universality of the ratio 
of entropy density and shear viscosity holds. This is because the shear viscosity depends on 
gravitational perturbations in tensor modes of SO (3) in the dual gravity. Each tensor mode 
satisfies a massless scalar field equation decoupled from the others, which ensures universality. 
The anisotropic symmetry-broken phase provides an anisotropic shear viscosity to entropy 
density ratio. The reason is that once 50(3) symmetry is broken to 50(2), the gravitational 
wave modes in the broken-symmetry directions are no longer 50(2) tensor modes. They are 
not decoupled from other fields and in fact interact with Yang-Mills fields. Therefore, in this 
case the gravitational modes do not display universal behavior. 
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Near the critical point (fx = fi c ), the phase transition is expected to depend on a 2 = %[14j. 
where k§ is 5-dimensional gravity constant and g is Yang-Mills coupling. If a is less than a 
certain critical value a cr u, the phase transition becomes second order. Near the second order 
phase transition, the ratio displays a scaling behavior with some critical exponent (3, 



where the value (3 = 1.00 ± 0.03 has been calculated numerically, taking into account the 
back reaction to the background metric. For a > a cr u, the phase transition is first order. 

3 This universality can be violated when string effects or quantum effects are to be taken into account in 
dual gravityPl9lfT0]. 

4 See [21 EQl [2TJ [22] for pioneering works on connection between Einstein-S'C/(2)-Yang-Mills and p-wave 
holographic superfhhds. 




(1) 



2 



Unfortunately, the large-a region has not been explored very well numerically, due to technical 
difficulties. Another interesting remark from [13] is that the critical exponent, (3 does not seem 
to depend on a. 

To examine such properties rigorously an analytic approach is needed. Our starting point 
for such an approach is a zero-mode solution at zero gravitational coupling, which is known 
exactly at the critical point [12]. In this note, we perturbatively analyze the properties of 
the anisotropic fluid near the critical point. We obtain the back reaction to the space time 
metric and also solve linearized equations of the gravitational perturbations and Yang-Mills 
fields from the back-reacted metric, using a double expansion of inverse Yang-Mills coupling a 2 
and eD%. e is a dimensionless small parameter and D\ is an 50(3) symmetry-breaking scale 
appearing in the anisotropic part of the Yang-Mills field. Similar perturbative expansions for 
different holographic models have been discussed by various authors [151 HI HEl [19]. However 
unlike those works, we consider both an anisotropy and the gravitational back reaction from 
the Yang-Mills fields. 

In our double expansion, the nontrivial leading order turns out to be 0(a 2 e 2 ). We also get 
the shear viscosity and entropy density ratio up to this order. Nonuniversality of the ratio of 
shear viscosity to entropy density shows up in the directions of broken symmetry while the 
ratio in the unbroken direction displays the expected universality. In our perturbative analysis 
we do not see the first order phase transition because a is small in the perturbative regime. 
Near the critical point, our solution also presents the scaling behavior as Eq([T]) and it turns out 
that the critical exponent f3 — 1 up to the leading order correction, which is consistent with 
the numerical result in [14]. In general, our perturbative results agree with and complement 
the numerical results in |14j . 



2 Holographic Setup and Large Coupling Expansion 

We consider the Einstein-5f/(2) Yang-Mills system in asymptotically AdS^, spacetime. The 
action is 

r / 1 1 o i 

? aMN 



/ 1 19 1 

d 5 xV^I^(R+ T2 )-— 2 F° MN F° 



(2) 



where M, N... are 5-dimensional space-time indices, a., are SU(2) indices and g is the Yang- 
Mills coupling. We conventionally choose L = 1. The Yang-Mills field strength F MN is given 
by 

Fmn = &MA a N — 9nA m — e a c A M A c Nl (3) 

where e abc is anti-symmetric tensor with e 123 = 1. The equations of motion from the action are 
obtained as 

Wmn = Rmn + 4Gmtv — K \ ^Tmn — -TpGMN^j = 0, (4) 
Y aN = V M F aMN - e abc A b M F cMN = 0, (5) 
where Tmn is energy-momentum tensor, of which form is 

Tmn = (^F MP F^ a — -F P Q a F PQa GMN S j ■ (6) 
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Our ansatz for the metric and Yang-Mills field are given by 

A = <p(r)r 3 dt + u(r)r 1 dx, (7) 
ds 2 = -N(r)a 2 (r)dt 2 + + r 2 f~ (r)dx 2 + r 2 f 2 (r) (dy 2 + dz 2 ) , 

where r a = y and s a are Pauli-matrices. The Yang-Mills field equations of motion in terms 
of the above ansatz are 

r X _ „ M i n | Ar) N'(r) J'(r)\ 4>\r)uj{r) 

^ ^( r ) + (?-f^)^( r )-ZWW^)=0 > 
\r cr(r) / r z N{r) 

and the Einstein equations are 



JV 2 (r)<r 2 (r) <r(r) r a{r) N(r) r N(r) N(r) a{r) N{r) 

2k 2 (f\r)u' 2 {r) + 2c// 2 (r) + 2^^' = (j 



3g 2 V. r 2 iV(r)a 2 (r) r 2 a 2 (r)7V 2 (r) 

2r 

^ = 2 . Vat? , W tt + 2r 2 N(r)W rr (10) 

J z (r) a(r) g z r \ N 2 {r)a 2 {r) J 

W = 2W yy + f(r)W xx (11) 

4r 2 N'(r) cr'(r) r 2 tS6 l2 (r) 

— 9 I- r — -I- r — -I = D 

N(r) ^ N(r) a(r) 3g 2 a 2 (r)N(r) ' 

f\r)W xx _ f"(r) | |V(r) | 3 | N'(r)\ f'(r) 1 | 2 JV'(r) 



2r 2 iV(r) /(r) V.^( r ) r N(r) J f(r) r 2 N(r) 2rN(r) 2ra(rj 

k 2 ( u' 2 {r)f\r) <p 2 (r)u 2 (r)f\r) c// 2 (r) \ J' 2 (r) _ 

^ ^ 3r 2 3rV 2 (r)iV 2 (r) 6iV(r)a 2 (r) / / 2 (r) 

W yy = W zz . (13) 

A known exact solution of the equations of motion is the AdS charged-black-brane solution 
given by 

2 

0(r) = Ml-^f), w(r) = 0, (14) 
<r(r) = /(r) = l andiV(r) = iV (r)=r 2 -^ + ^^l, 

where /z is chemical potential, is the black brane horizon and m = H ^ — -. In the 

infinite Yang-Mills coupling limit as g — > oo, the last term in iV(r) vanishes and the solution 
becomes uncharged. 
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2.1 Large Coupling Expansion and its Leading Order Corrections 

In this section, we develop corrections to the metric and Yang-Mills field perturbatively in a 
double expansion in sD\ and a 2 = . £ is a dimensionless small parameter and D\ is the 
(SO (3) rotational symmetry-breaking order parameter. We choose the horizon of the black 
brane to be conventionally located at r = 1 by scaling r — > r^r and {t, x, y, z} — > ^r{t, x, y, z} 

and defining a new chemical potential // = The equations of motion enjoy a certain 
scaling symmetry [T3| [IT] . By means of these rescalings, we can choose the asymptotic values 
of oir = oo) = 1 and f(r = oo) = 1 at the large r boundary where the spacetime becomes 
asymptotically AdS§. The value of the chemical potential in the dual boundary field theory is 
taken to be /z = 4 at the phase transition point. To obtain corrections, we expand any fields 
a(r) appearing in the ansatz ([7j) as 

a(r) = ao(r) + eai(r) + e 2 a2(r)... (15) 

Each term in the expression can in turn be expanded as 

di(r) = a ij0 (r) + a 2 a i)2 + a 4 a iA (r)... (16) 

The zeroth-order solution in e is given in Eq ([T4"|) . where only N contains a subleading correction 
of order a 2 in the sense of the above expansion. N Qj2 = y (pr — ^) and the higher-order terms 
in a 2 vanish, N ^ = for i = 4, 6.... The detailed computations of the nontrivial leading order 
corrections to metric and Yang-Mills field are given in Appendix [AT Here, we briefly list the 
leading-order back-reaction corrections to the metric, which are given by 

- ^so^p /w ^- £V lw (17) 

and w(r) = r > * + ^! C 1 _ 1 )_ eW ^l ( 1 + 2r* 



r 2 3 yr 4 r 2 J 9r 2 \r 2 (l + r 2 ) 3 2(1 + r 2 ) 2 

281 / 1_ 

+ 560 V r 2 

Any subleading corrections to the Yang-Mills field in a 2 would not contribute to the leading 
back-reaction corrections to the metric (Our aim is to get metric corrections up to 0(a 2 e 2 )). 
Therefore, we obtain the Yang- Mills field solutions up to fa o and Ui o only. These are given by 



"(0 = e 7^3 + 0(^ 



1, e 2 D 2 ( (l + 2r 2 ) 1 281 / 1 \\ n . „ 
The black brane temperature is modified by the leading corrections to 

r^-T^m*™)' (20) 
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where T c = - (l — |a 2 ) is the critical temperature at the phase transition from the isotropic 
phase to the anisotropic phase. The black brane entropy is 

S = ^ (21) 

where V3 is spatial coordinate volume of the boundary space-time, V3 = J dxdydz, in this 
rescaled coordinate. 



3 Anisotropy of Shear Viscosities 

In this section, we calculate the ratio of shear viscosity to entropy density via the Kubo formula, 
by considering fluctuations Hmn and 5A a M around the background metric and the background 
Yang-Mills field, respectively. We choose the gauge hur = $A a r = 0. In the anisotropic 
phase, the bulk gravity system enjoys residual SO(2) and Z2 symmetries. The modes may be 
decomposed according to their SO{2) representations as 

• Scalar modes in SO{2) : h yz ; h yy — h zz , 

• Vector modes in SO{2) : h yt ,5A 3 y , h xy ,8Al,5A 2 y ; h zt , SA 3 Z ; h xz ,8A\,5A% 

• Tensor Modes in SO{2) : h tt ,h yy + h zz ,h xx ,h xt ,SAl, 5A a £ , 

where each decoupled mode is categorized by semicolons. h yz is totally decoupled from any 
other modes and satisfies a massless scalar field equation showing universality. However h xy 
interacts with 5A y and 5 Ay, leading to nonuniversal behavior. In the following, we will obtain 
solutions for h yz and h xy and show this explicitly. Other modes can be calculated by the similar 
methods. 



3.1 Universality of 

In this subsection, we calculate the ratio of shear viscosity to entropy density for the shear 
mode h yz , using the double expansion that we introduced in the previous section. We consider 
fluctuations of the Yang-Mills field and metric fields around the background metric (TTTT) and 
obtain perturbative corrections up to 0(e 2 a 2 ). We only consider time dependent fluctuations 
with frequency v and use small frequency expansion up to first subleading order in v. Even 
in the presence of a nonzero symmetry breaking parameter JDj, the rotational symmetry in 
y — z plane is not broken and the ratio — is universal. In the following, as a warm up, we will 
explicitly calculate this ratio to be j- up to leading-order corrections in a 2 e 2 . To show this, 
we begin with the linearized equation of motion of h yz {r) = r 2 / 2 (r)$(r, t), 

= *» + (- + -^r - ^^k) 4» + (22) 



r N(r) 3(7 2 (r)N{r)J vK ' N 2 (r)a 2 (r)' 

where the prime denotes the radial derivative. For the field $(r), we have used the Fourier 
transform from real time to frequency, as 

POD 

$(r,t) = / e~ iv % v {r)dv (23) 
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The near horizon behavior of <& u (r) should be a purely ingoing solution 



r 



1\ -*l (*+¥ ^-^o^-K^o*) 

$„(r) ~ ( 1 - - , (24) 



where A; and / are integers with k > 2 or I > 2. With this boundary condition, the solution 
$y(r) is obtained as 



iV(r) 



-^(l + f Q 2 -^ £ W?)+0( £ V) 



^(0=(^) ^,« 2 ), (25) 

where 



oo 



F(e,a 2 )= ^$, 2j (r)£^. (26) 

i,j=0 

Each $j ( 2j and its near-AdS" boundary expansion are given in Appendix [Bl Here, we briefly 
discuss the near boundary expansion of the solution to get — . Defining the boundary value 
$ = $„(oo), $j/(r) can be expanded as 

$ v (r ->■ oo) = $ + + 0(rVe*a') (27) 

in the large r limit (See Eq (l7ip in Appendix IB"]) , where i < —4, j > l,k > 2 or I > 2. Using 
the prescription to get retarded green's function in [T4], we get 

G^ yz (u,k = 0) = - Y + O(u 2 ) (28) 



2/cg 



The shear viscosity in y — z direction is given by 



^^h [G "'- G ^ = ^ (29) 

where star indicates complex conjugate. Using the entropy of the black brane( l2"TT) . the ratio of 
shear viscosity to entropy density is obtained as 

^ = — . (30) 

S 47T 



This value turns out to be universal up to 0{e or). 



3.2 Nonuniversality of ^ 



We start with a set of equations with h xy = r 2 f 2 (r)^/(r,t), 5A y and 5 A 2 . The superscripts 
on 5 A fields note SU{2) indices and subscripts do space-time indices. We also solve these 
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equations in the frequency space by Fourier transform as in the previous subsection. The 
equations in the frequency space are given by 

„, (1 _jb_ Qf'jr) ra 2 (j)' 2 (r) \ u 2 ^(r) 
' ~ \r + N(r) + f(r) 3N(r)a 2 (r)J W + N 2 (r)a 2 (r) { °- } 

2a 2 f , v _ u(r)<p 2 (r)5Al(r) tuu(r)<j)(r)5A 2 y (r) 
r 2 f 2 (r) V [r) y[T) N 2 (r)a 2 (r) N 2 (r)a 2 (r) 



+ 



r f(r) N(r) a(r) J yy/ \N 2 (r)a 2 (r 
f(r)u\r)^\r) - 



2iis(f)(r)5A 2 y (r) 



N 2 (r)a 2 (r) ' 



r f(r) N(r) a(r) J v \N 2 (r)a 2 (r 



We expand each field with the same fashion as Eq(f25 



N(r)\ °?-&5**<*i%)+o&« 1 ) 



*(r) = GCe,^), (34) 



.U, L ,(/-) = (^] HC.d 



3 

, J, 



_ i |(l + ^ a 2__ll_ £ 2 a 2o2) +0 ( e fc a i) 

2 ~ 



where the functions G(e,a 2 ), H(e,a 2 ) and I(e,a 2 ) are expanded as 

oo 

G(e,a 2 ) = ^TVwWe'aF, (35) 

i,j=0 
oo 

H(e,a 2 ) = 8A^(r)JaF, 



ij=o 



OC 



and I(e,a 2 ) = V] <L4? 2y (r)g*Q! 2j '. 
i,i=o 

The solutions of the equations are listed in AppendixO Corrections to the Yang-Mills 
fields which are subleadings in a 2 do not contribute to the leading order corrections of the 
shear viscosity, so we get 5Aj (r) and 5A 2 (r) only(See Eq (l3~Tj) ). We also use small frequency 
expansion as in the last subsection, and obtain the solutions up to 0{y). We specify purely 
ingoing boundary conditions for the fields, the form of which are the same as Eq (124p . Here, 
we discuss the near boundary expansion of *&(r), obtaining the retarded correlator of h xy and 
the shear viscosity to entropy density ratio, — . 



S 



The near- AdS boundary expansions of ^(r), 5A y and 5 Ay are given by 

Mj(r = oo) = ~(6A^-22A^-D 1 ^)+0{e^) i (36) 
<L4j(r = oo) = ^(22Af} ) -QA^-llD l ^) + 0(s,% 



and 



*M = (42+^S + £ 2 ^)+a 2 (4?+ e <2+^) (37) 

The .SO (3) symmetry is broken spontaneously, so the Yang-Mills field should not provide any 
source terms to the dual field theory system. Therefore, A y (r) and Ay(r) should become 
normalizable modes of the solutions, then we have 

6ig - 22ig> - Dnf,® = 0, and 22A% - QA§ - UD^ = 0. (38) 

The solutions of these equations are 

A§ = and A% = (39) 

Using Eq.(l39l) and defining \l/(oo) e f as a boundary value of *&(r), the near boundary 
expansion of \l/(r) is given by 

tt(r -> 00) = v& + — vl/ + eW^^. (40) 

The prescription of the retarded green's function in [T3] provides 

-%v 29 
T2 

'5 

and the shear viscosity is calculated as 



G%Jy, k = 0) = ^ I 1 + g^VC? ) + 0(„ 2 ) ( 1 1 ) 



*• s as 2is' G ™ - G -»> = 2^i( i+ w/ a2i)2 ) ■ (42 > 

Using entropy of the black brane (l2~Tj) . the ratio of shear viscosity and entropy density obtained 



as 



Vxy _ 1 / , . 29 o 2 n- 



Therefore, the shear viscosity and entropy ratio in x — y direction is not universal, and we have 
shown this up to non-trivial leading order correction in a and e. Using the temperature of the 
black brane, Eq (|43p can be written as 



j] xy 1305ttT c ( a T 



= _____ _ 

s 544 V T cJ 

where /3 = 1. It is also shown that the critical exponent j3 = 1 up to corrections of order e 2 a 2 , 
near the phase transition point T = T c . 



Appendix 

A Leading Order Solutions 

In this section, we solve EqfllJ) and Eq(J5]) to get leading order back reaction to the metric 
by symmetry breaking order parameter in Yang-Mills field, D\. Using metric and Yang-Mills 
fields ansatz© and their expansions ( !T5|) and (Fl6|) . we obtain the equations order by order in 
e and a. Since the aim is to get leading back reaction to the metric, we get Yang-Mills field 
solutions up to the zeroth order in a and the first order in e for uj(r) but the second order in 
e for <p(r), both of which provide first leading order corrections to the metric correctly. The 
first order Yang-Mills equations in e are given by 

= #(r) + -#(r)- Vi( r )> ( 45 ) 

... . (\ 2(r 4 + 1) +r 3 N'(r)\ ... l&^fr) , s 

= <4(r) + - + „ , ° _ u[(r) + — - )> — (46) 



r r(r 4 - 1) +r 3 N (r)J 1W ( r 2 + i + ^y) )2 
and the first order Einstein equations in e are 

= oKr), (47) 

64a 2 /r 3 \ 
= 2 ^(r)+r^(r) + -g£- (_#( r ) - 2^(7-)) , (48) 

= 2rN 1 (r) + r 2 N[(r) + 2(1 -5r 4 )/((r) + 2r(l - r 4 )/f (r) (49) 
16a 2 

(-2r 2 (2JVi(r) + WV((r)) - 8r(l + r 2 - 2r 4 )/{(r) - r 3 (l + r 2 )4>[(r)) . 



3r 3 (l + r 2 ; 



We expand every field appearing in above equations using Eq (fl6|) . With this expansion, 
the solutions of metric are trivial, which are given by 

ci,o(0 = ^1,0, <ri,2( r ) = °> ( 50 ) 
/i,o(r) = /i,o, /i, 2 (r) = 0, 

^,o(r) = % N h2 (r) = 8 -cJ 1 ^ 
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where we set a^o = fi t o = for the boundary values cr(oo) = f(oo) = 1. N\ t o = and = 1 
for the space-time has its horizon at r = 1. The Yang-Mills fields equations (I46p up to 0(a 2 ) 
are given by 

n " f \ i 1 + 3r4 / f \ , 16w i,oW / -i \ 

= M r ) + r(r 4_ 1 ) a; l,o( r )+ ( r 2 +1 )2 ' ( 51 ) 

n // / \ , l + 3r 4 , lft^ 128A(l + 9r 2 -2r 4 ) 

= u (r) + — — -(J u (r) + + 



r(r 4 -l) Mvy (r 2 + l) 2 3r 2 (l + r 2 ) 5 

The Yang-Mills field solutions are solved as 

D x r 2 

= ( r 2 + l)2 ' ( 52 ) 



^i,o(r) = yf 1 -^)' (53) 



P 2 r 2 128Ar 2 f (y 2 + l) 3 / 136 1, ( 2y 2 , 
7 + 53y 2 + 43t/ 4 + 27y 6 + 6y 8 



12(1 + y 2 ) 5 

The second order equation of Yang-Mills field in e is 



3 8 4r 2 D 2 

= W + -0' 2iO (r) - -4 (r) - (55) 



and Einstein equations are obtained as 



. , . { o 1 32a 2 32ct 2 \ , . 

= 2N 2 (r)+rN>(r) + 2a>(r)(r 2 --- — + —) (56) 

^fM + T 0' 2 (r)J , 



64a 2 /V 

+ 



3r 4 

, . , 8a 2 D 2 r , , 

= 6a' (r) -— -, 57 

2V ) 3 ( r 2 + 1 )4' v ; 

aw,„ x , *w„ n 4a 2 D 2 r(r 2 -l)(l-6r 2 + r 4 ) , . 

= r(-l + r 4 )/^(r) + (-1 + 5r 4 )/^(r) — ^ . \ ... (58) 



(1 + r 2 ) 5 

where the second order equation of ui(r) is not given, because which provides sub leading 
corrections to the metric backreaction. 

The Yang-Mills field solution of (p(r) in 0(e 2 ) are given by 

= 02,o + C 2 (l - + ^(1 + ^ (59) 
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and the metric corrections are 



O"2,0 


= &2,0, 




2-L^ / 

9 V 


N 2 , (r) 


N 2 ,o 


N 2 , 2 (r) 


16 | 


Mr) 


— /2,0, 


/a >2 (r) 


- f ' 

J2,2 



(60) 



(1 +r 2 




18(1 + r 



2^|4 



(62) 

2 / l+2r 2 3r 2 \\ 

r 2 12 V^ 2 (l + r 2 ) 3 2(l + r 2 ) 2 ; y /' 1 J 

(64) 
(65) 



D 2 ~ 

where 02,o = — "32" f° r the regularity of the Yang-Mills field at the horizon and N 2y0 , f 2 ,o, ^2,2 
and f 2}2 are 0(1) constants which are set to be vanished for /(oo) = cr(oo) = 1. iV 2i 2 =32+02 
for the space-time has its horizon at r = 1. C\ and O2 are coefficients of zero modes of Yang- 
Mills field equations. Without loss of generality, we can set C\ = 0. However, C 2 = 1 
requesting u)2,o(r) to be normalizable mode and regular at the black brane horizon. Then, we 
have metric backreaction only with D\ as a SU (2) symmetry breaking scale. 

B Leading Order Correction of h yz 

In this section, we briefly describe the solutions of Eq (l22]) using the form of expansion ( 1251) and 
(1261) . Each term in expansion ( 1261) is given by 



*o,o(r) = C + ^S + O^ 2 ), (66) 
*o, 2 (r) = ^l + ^(^l + ^ol(ln(l + \) -\)) +0(u% (67) 



^2 / ^->2 

$i,o(r) = 0S + ^ + O(z/ 2 ), $ 2iO (r) = 0^ + z/0« + O(z, 2 ), (6 



$ li2 (r) = 0S + ^(0S + 8#S(^(l + - --))+0(z, 2 ), (69) 



and 



,(0) . ,(i) if / ,(0) f, 2 / 105 175 70 



* 2 , 2 (r) = C + ^-^^x 2 ^^ (70) 
35 279\ , 67200$ ino ^(o) «.<D) 



192 D{£ -35^ Zn(l 



(1+r 2 ) 4 2r 2 7 r 2 ^—j.tu.u ^y 2 ,uj -■"\^ 1 r2 , 

We also obtain the near AdS boundary expansion of $(r) as 

$,(r^oo) = <^e$a + , 2 $g + a 2 ($g + ^S + ^S) (71) 

12 



+ 0(r i u j e k a l ), 
where z < —4, j > 1, > 2 or I > 2. 

C Leading Order correction of h xy 

In this section, we list the solutions of the set of equations AST]) . (152"]) and (|33|) . We listed our 
solution using the expansion (134")) and ( 135|) . As explained in Sec J3.2[ we get 8A} and <L4? only 
for the Yang-Mills field solution. 5Aq and 5A^ are zero modes of the solutions. Without 
loss of any generality, we set 8A\ = 8A% = 0. The first subleading corrections of Yang-Mills 
fields in e are given by 

KoW = (T^^S + ^o- J^{^(l + r* + 24rHn(r) (72) 

- 8r 2 /n(l + r 2 )) + 2^(5 - llr 4 - 40r 2 /n(r) - 24r 2 /n(l + r 2 )) 

- L>i4°o(17 + r 4 + 56r 2 /n(r) - 24r 2 /n(l + r 2 ))) + 0(z/ 2 ) 

and 

KoW = (^°o + + l|^H o(- 5 + U ^ + 40r 2 Mr) (73) 

+ 24r 2 /n(l + r 2 )) - 6^(1 + r 4 + 24r 2 /n(r) - 8r 2 /n(l + r 2 )) 

+ Di^ ] (5 - llr 4 - 40r 2 /n(r) - 24r 2 /n(l + r 2 ))) + 0(z/ 2 )) . 

\I/(r) solution is also obtained as 

*o,o(r) = + ^ + (74) 
*i,o(r) = ^3 + W$ + C>(^), (75) 
* 2l o(r) = ^°,] + ^S + 0^ 2 ), (76) 

*o, 2 (r) = ^ + ^+8^^n^l + ^-^+0(^), (77) 

and 



*i, 2 (r) = Ci + ^ + Sz^ (ln(l + -) --) + 0(z/), (7E 



W = « + (if^-iiif^ (79) 
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+ 



576 \ (1 + r 2 ) 3 

2«£ 1 (3Ag> - 5jgS + 4gijg) 2iD 1 (-9i$ + 19j^ + 37giVg) 
1 + r 2 (1 + r 2 ) 2 

12D!(3ii^ + 48A$ - 3z4°i - ^i4°o) 36i(93D 2 4°o - 4480^$) 



(1 + r 2 ) 4 



(1 + r 2 ) 4 35r 2 

18ig + 10A® + (80) 



8iD x (2r 2 - l)/n(r 
-5A 

)H 



- 4i(-3AS3A + SlgiA - + 2304^)/n(r) 

35 



+ tf„(l + r 2 )(-6iS + 10^2 - ^i^Vffl + 460 <o) 
+ 0(^ 2 ). 
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